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ABSTRACT 


The  C-vector  method  is  used  to  determine  the  three-dimensional 
pseudo- vorticity  field  of  the  East  Greenland  Current  by  using  a  CTD 
data  set  acquired  during  the  1984  Northwind  cruise.  The  value  of 
the  ageostrophic  pseudo- vorticity ,  the  total  (geostrophic  and 
ageostrophic)  psaudo-vorticity  and  the  vertical  component  of  the  C- 
vector  curl,  can  be  calculated  by  the  C-vector  method.  From 
these  values  the  pseudo- vorticity  can  be  investigated  for  the  East 
Greenland  Current.  Additionally,  the  positive  and  the  negative 
vertical  motion  at  each  level  can  be  investigated  from  the  'ir  value. 
A  two- cell  structure  with  downward  motion  around  the  East  Greenland 
Polar  Front  was  found.  The  C-vector  method  also  demonstrates  that 
the  vertical  circulation  induced  by  an  anticyclonic  gyre  and  the 
cross -coastal  circulation  generated  by  the  surface  wind  can  be 
depicted. 
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I. 


INTRODUCTION 


A.  GENERALITY 

It  has  been  a  long-term  effort  for  physical  oceanographers 
to  determine  the  three  dimensional  (3-D)  circulation  of  the 
ocean  based  upon  measurements  of  its  internal  density 
distribution.  Two  distinct  approaches  have  been  developed; 
(1)  the  dynamical  method,  which  is  widely  used  by 
experimentalists,  in  which  one  computes  the  dynamical  height 
of  the  sea  surface  relative  to  a  certain  depth,  from  which 
geostrophic  currents  can  be  calculated  by  assuming  a  certain 
level  of  no  motion;  and  (2)  the  numerical  method,  mostly  used 
by  modelers,  which  solves  the  complete  prognostic  dynamical 
equations  under  some  assumptions  such  that  the  density  field 
does  not  change  with  time.  The  dynamical  method  computes  only 
the  geostrophic  current;  no  3-D  circulation  can  be  inferred  by 
this  method.  However,  the  numerical  method  of  computing  the  3- 
D  ocean  circulation  based  upon  an  observed  hydrographic  data 
set  presents  a  difficult  problem,  i.  e.,  the  numerical  method 
does  not  provide  a  one-to-one  correspondence  of  the  3-D 
circulation  to  the  observed  hydrographic  data  set.  In  order 
to  overcome  this  difficulty,  the  C-vector  method  is  proposed 
in  this  thesis. 
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B.  TRADITIONAL  METHODS  OF  CURRENT  CALCULATION 


Traditionally,  two  methods  exist  co  investigate  oceanic 
currents,  the  direct  method  and  indirect  method.  Before  the 
early  years  of  this  century,  scientists  put  drift  bottles  in 
the  ocean.  Some  of  these  bottles  were  found  and  returned. 
This  is  the  simplest  direct  method;  however,  it  provides 
trajectories  of  only  the  surface  currents.  Another  direct 
method  is  to  moor  current  meters  at  fixed  positions.  This 
method  can  obtain  current  data  for  several  different  depths. 
However,  this  data  pertains  only  to  the  immediate  vicinity  of 
the  current  meter. 

The  indirect  method,  also  known  the  geostrophic  or  dynamic 
method,  requires  information  of  the  density  distribution  in 


the  ocean. 

It  is  easier 

to  obtain 

this  information 

from 

measurements 

of  temperature  and 

salinity 

than 

from 

measurements 

of  currents 

directly. 

In  fact. 

most  of 

our 

knowledge  of  the  ocean  circulation  below  the  surface  has  been 
obtained  in  this  way. 

The  geostrophic  method  uses  the  density  field  in  the 
interior  of  the  ocean  to  calculate  the  horizontal  current 
field  by  the  geostrophic  equation.  If  the  slopes  of  constant 
density  (or  isopycnal  surfaces)  vary  with  depth,  this 
indicates  that  the  currents  vary  with  depth.  We  say  that,  in 
this  case,  there  is  'current  shear'  or  'geostrophic  shear’. 
This  is  characteristic  of  the  baroclinic  situation  in  which 
the  slopes  of  the  isopycnals  vary  because  the  water  density 
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depends  on  water  properties  {temperature  and  salinity)  as  well 
as  on  pressure  (depth)  .  The  barotropic  situation  is  one  in 
which  the  density  depends  on  depth  only;  there  will  be  no 
current  shear,  i.e.,  the  current  will  be  the  same  from  top  to 
bottom.  This  is  the  barotropic  part  of  geostrophic  current  and 
the  geostrophic  equation  gives  us  no  information  about  it. 
The  barotropic  geostrophic  currents  must  be  inferred  from  the 
trace  method  (Reid,  1989)  or  a  reference  level  (level  of  no 
motion) .  The  dynamical  method  only  resolves  the  geostrophic 
currents;  no  information  about  the  vertical  circulation  can  be 
obtained  from  such  methods. 

C.  THE  PURPOSE  OF  THIS  THESIS 

In  this  thesis  the  author  employs  a  method  of  calculating 
the  C-Vector  that  is  developed  from  the  geostrophic  method  in 
order  to  obtain  more  information.  The  general  idea  is  set  out 
in  meteorology  by  the  Q-vector  concept  (Hoskins  et  al. ,  1978)  . 
Hoskins  et  al .  determined  that  the  Q-vector  derived  from  the 
CO  equation  [the  (co)  equation  is  a  very  useful  diagnostic 
ecjuation  used  in  mid- latitude  cyclogenesis  studies  as  the 
vertical  velocity  is  important  for  understanding  the  observed 
cloud  structure  and  the  development  of  cyclones  and  fronts.  It 
describes  the  vertical  velocity  distribution  that  results  from 
dynamical  forcing.]  can  be  explained  on  the  synoptician's 
chart  more  precisely  than  earlier  studies  (Sutcliffe,  1947) . 
This  theory  is  very  helpful  for  learning  about  and  diagnosing 
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the  synoptic  and  frontal  vertical  circulation.  A  disadvantage 
of  this  concept  is  that  it  ignores  the  barotropic  part  of  the 
rotational  ageostrophic  flow.  According  to  Xu  (1992) ,  this 
problem  can  be  overcome  by  a  new  concept,  the  C- vector 
concept,  which  provides  a  method  of  resolving  the  ageostrophic 
circulation  conveniently  and  a  new  method  of  solution  which 
relies  on  the  C-vector  equation.  The  C-vector  concept  also  can 
be  applied  in  physical  oceanography  (Chu,  1992) .  This  thesis 
employs  the  hydrographic  data  that  was  collected  from  the  East 
Greenland  Current  by  the  icebreaker  ''Northwind''  in  1984.  A 
three-dimensional  model  will  be  designed  to  distribute  the 
irregularly  spaced  hydrographic  data  on  each  grid  point.  The 
C- Vector  will  then  be  computed  from  these  gridded  data. 

D.  ENVIRONMENTAL  CONDITIONS  IN  THE  EAST  GREENLAND  CURRENT 
1.  Water  Hasses 

Four  primary  water  masses  have  been  recognized 
historically  within  the  East  Greenland  Current  (EGO  north  of 
the  Denmark  Strait  and  in  the  Greenland  Sea  (Coachman  and 
Aagaard,  1968  a,b) .  They  are  known  as  Polar  Water,  Atlantic 
Water,  Atlantic  Intermediate  Water  and  Deep  Water.  Although 
these  water  masses  are  well  recognized  features  of  the 
Greenland  Sea,  their  properties  have  been  variously  defined  by 
polar  researches.  The  following  discussion  of  water  mass 
properties  follows  that  of  Paquette  et  al.,  1985  and  Bourke  et 
al.,  1987. 
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a.  Polar  Water 


Polar  Water  (PW)  is  defined  at>  having  a  temperature 
of  less  than  0°C  and  a  salinity  of  less  than  34.5  psu.  Water 
with  these  characteristics  occupies  the  upper  layer  of  the 
EGC,  and  extends  over  the  east  Greenland  shelf,  reaching 
depths  of  over  200  m  in  places  (Figure  1.1).  Near  the  East 
Greenland  Polar  Front  (EGPF)  it  is  pinched  off  at  the  surface 
by  the  warm  water  of  the  Return  Atlantic  Current (RAC)  to  the 
east . 

h.  Atlantic  Water 

Atlantic  Water  (AW)  has  its  origin  in  the  warm, 
saline  Norwegian  Atlantic  Current  and  has  salinities  above 
34.9  psu  and  temperatures  above  3°C,  sometimes  exceeding  5°C 
even  at  latitudes  north  of  75 ®N. 

c.  Atlantic  Intermediate  Water 

Atlantic  Intermediate  Water  (AIW)  ranges  in 
temperature  from  0®C  to  3°C  and  in  salinity  from  34.5  to  34.9 
psu.  This  water  is  found  over  the  deeper  parts  of  the  shelf 
and  the  troughs  which  cut  across  it  and  in  the  region  of  the 
EGPF.  It  is  also  found  in  the  upper  layer  between  the  West 
Spitsbergen  Current  (WSC)  and  the  EGC.  In  the  frontal  zone  and 
westward  thereof,  it  lies  underneath  the  PW.  The  RAC  is 
essentially  all  AIW.  The  RAC  has  its  source  in  the  AW  of  the 
WSC,  part  of  which  turns  westward  and  southward  at  latitudes 
between  75°N  and  80°N  or  occasionally  even  farther  north. 
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d.  Deep  Water 

Deep  Water  (DW)  is  found  below  about  300  m.  The 
temperature  of  this  water  is  less  than  0°C  decreasing  slowly 
until  it  attains  a  value  near  -i°C  at  the  bottom.  The  salinity 
varies  little  remaining  between  34.88  and  34.90  psu,  which 
represents  a  decrease  from  the  usually  higher  salinities  of 
the  AIW  above  it.  No  DW  is  found  on  the  shelf. 

2.  Currents  and  Circulation 

a.  General  Circulation  in  the  Greenland  Sea 

The  general  circulation  in  the  Greenland  Sea  is 
presented  in  Figure  l.l.  The  surface  circulation  is  dominated 
by  a  large  cyclonic  gyre  bounded  to  the  south  by  the  Jan  Mayen 
Current  and  to  the  east  by  the  Norwegian  Atlantic  and  West 
Spitsbergen  Currents.  In  the  northern  portion  of  the  Greenland 
Sea,  the  West  Spitsbergen  Current  splits  into  two  branches, 
one  branch  continues  to  flow  northward;  the  second  turns 
westward,  subsequently  submerging  to  about  200  m  and  turning 
southward.  This  second  branch  of  relatively  warm  water,  called 
the  Return  Atlantic  Current  (RAC) ,  together  with  the  East 
Greenland  Current  (EGC) ,  which  comes  from  Arctic  PW,  accounts 
for  the  flow  of  water  in  the  western  portion  of  the  Greenland 
Sea.  Figure  1.2  is  a  map  of  the  estimated  near  surface 

circulation  pattern  in  the  region  of  the  East  Greenland 
Current  based  upon  the  hydrographic  data  obtain  during  the 
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Figure  1.1  General  circulation  in  the  Greenland  Sea  (from 
Paquette  et  al . ,  1985). 

1984  Northwind  cruise.  During  the  period  of  this  cruise  two 
current  meters  were  moored  on  the  1000  m  isobath  near  78.5°N 
essentially  under  the  jet  of  the  EGC  (Muench  et  al . ,  1986). 
The  mean  speed  compared  well  with  the  baroclinic  current 
calculated  by  Bourlce  et  al.  (1987),  based  on  geostrophic 


Figure  1.2  Near-surface  baroclinic  circulation  over  the 
east  Greenland  continental  shelf.  The  length  of  the  arrows 
is  suggestive  of  the  relative  speed,  with  maximum  speeds 
near  the  EGPF  (Bourke  et  al.,  1987) 
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velocity  cross  sections  referenced  to  500  m  depth  or  the 
bottom.  The  southward  part  of  the  geostrophic  Transport 
ranged  from  1.33  to  1.61  Sv,  with  an  average  of  1.47  Sv.  From 
three  sections  which  extended  across  the  entire  shelf,  the 
northward  part  of  the  transport  ranged  from  0.44  to  0.72  Sv, 
with  a  mean  of  0.58  Sv.  These  means  were  subtracted  to  give  a 
net  flow  in  the  EGC  of  0.89  Sv  southward.  The  northward  flow 
near  the  coast  turned  anticyclonically  near  80°N  following  the 
depression  of  the  Westwind  Trough.  This  anticyclonic 
circulation  transported  about  0.5  Sv,  providing  a  substantial 
amount  of  recirculating  upper  layer  shelf  water  to  the  frontal 
zone.  The  northward  flow  near  the  coast  also  was 
substantiated  by  ice  drift  measurements  using  NOAA  7  imagery 
in  August  -  September  1985. 

b.  Return  Atlantic  Current 

The  Return  Atlantic  Current  (RAC)  is  found  to  be  a 
notable  feature  and  requires  more  emphasis  than  that  given  in 
previous  work.  The  current  gains  its  identity  from  being  a 
core  of  warm,  high- salinity  water,  often  broken  into  differing 
filaments,  most  of  it  being  above  O^C  and  near  the  high  end  of 
the  salinity  range  for  AIW,  even  overlapping  into  AW.  Property 
limits  for  these  waters  are  described  later  in  this  paper.  The 
RAC  has  a  width  of  100  km  or  less  and  lies  along  the  EGPF  at 
depths  roughly  between  50  to  300  m,  commonly  at  least  partly 
under  the  PW. 
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3 .  East  Greenland  Polar  Front 

This  front  acts  as  the  boundary  between  cold  fresh 
PW,  flowing  southward  as  part  of  the  EGC  along  the  Greenland 
coast,  and  the  warmer,  more  saline  AIW  immediately  to  the  east 
(Paquette  et  al.,l985)  .  The  warm,  saline  water  to  the  east  of 
the  ice  must  originate  in  the  AW  of  the  WSC.  At  the  ice  edge 
the  horizontal  density  gradient  (i.  e.,  the  baroclinic  shear) 
is  strongest.  It  is  in  this  region  that  a  narrow  (~50  km  wide) 
southward  flowing  jet  exists.  Maximum  speeds  in  excess  of  40 
cm/s  have  been  observed.  The  jet  decays  with  depth 
essentially  vanishing  by  300  m. 


break  (Bourke  et  al.,  1987). 
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Figure  1.3  is  a  temperature  salinity  cross  section 
along  78°12  N.  In  the  vicinity  of  station  199,  a  sharp 
temperature  and  salinity  gradient  can  be  seen.  This  feature 
marks  the  core  of  the  EGPF.  The  cold  and  dilute  PW  can  be  seen 
to  the  left  of  the  EGPF;  the  warm  and  saline  AW  and  AIW  are  to 
the  right  of  the  EGPF. 

Figure  1.4  is  a  similar  cross  section  farther 
southward  aligned  along  the  axis  of  Belgica  Trough.  In  this 
figure,  the  region  of  the  sharp  gradient  of  temperature  and 
salinity  is  between  stations  265  and  263. 

The  data  in  these  two  figures  represent  summertime 
conditions  prevalent  during  1984, 

E.  DATA  BACKGROUND 

1.  Original  Hydrographic  Survey  Data  and  Grid  Point  Area 
a.  Original  Hydrographic  Survey  Area 

During  August  and  September  1984,  the  icebreaker 
Northwind  visited  the  northern  part  of  the  Greenland  Sea,  the 
east  Greenland  shelf  and  the  EGPF  region.  A  hydrographic 
survey  was  conducted  in  the  region  bounded  by  IV^W  and  10°E, 
and  TB.S^N  and  81.5*N.  More  than  300  CTD  stations  were 
occupied  during  this  four  week  cruise  (Figure  1.5) .  This 
study  uses  the  data  from  23  7  stations  from  this  cruise  (Figure 
1.6).  In  addition  to  temperature  and  salinity  data,  density 
(Cj)  and  wind  speed  are  also  used. 
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axis  of  Belgica  Trough  and  across  the  EGPF  (Bourke  et  al .  , 
1987)  . 

b.  Grid  Point  Area 

The  region  encompassed  by  the  Northwind  cruise  was 
gridded  to  provide  a  spatially- smooth  data  field.  The  grid 
extended  from  IS^W  to  S^W  and  from  to  Sl^N.  The 

rectangular  grid  area  contains  10  x  21  grid  points  with  the 
distance  between  each  point  being  about  25  )an.  The  northwest 
corner  of  this  grid  area  is  on  land.  Excluding  the  points  on 
land,  the  grid  contains  195  points  (Figure  1.8). 

2 .  Bathymetry 

Johnson  and  Eckhoff  (1966)  collected  previously 
acquired  bathymetric  data  from  the  Greenland  Sea  to  make  a 
bathymetric  map  of  that  region.  According  to  their 
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Figure  1-6  The  stations  < 
cruise  for  distributing  the 


Figure  1.7  Wind  vectors  at  CH)  stations  acquired  during  the 
Northwind  1984  cruise. 
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Figure  1.8  Positions  of  grid  points  in  the  three- 
dimensional  C-vector  model  of  the  western  part  of  the 
Greenland  Sea. 

description,  the  Greenland  shelf,  as  is  typical  with  most  high 
latitude  continental  shelves,  is  broad,  irregular,  and  fairly 
deep.  The  continental  shelf  is  approximately  100  km  wide  at 
70°N  and  broadens  to  240  km  at  76®45’N.  It  then  narrows  with 
increased  latitude  being  only  150  km  wide  at  80°30’N.  The 
continental  shelf  is  typically  rough  with  an  approximate  9-m 


amplitude  and  700 -m  wave  length  relief  superimposed  on  a 
rather  uneven  platform.  Troughs  are  present  on  shelves  near 
the  coast.  Some  of  these  appear  to  parallel  the  coast  of 
Greenland  while  others  are  nearly  normal  to  the  shoreline. 

A  detailed  map  of  the  east  Greenland  shelf  bathymetry 
is  shown  in  Figure  1.9.  It  is  derived  from  the  map  of  Perry  et 
al .  (1980)  but  significantly  updated  by  depth  measurements 

made  at  CTD  stations  during  the  Westwind  1979  and  the 
Northwind  1981  and  1984  cruises  and  produced  by  Tunnicliffe 
(Bourke  et  al.,  1987).  The  bathymetry  of  the  Greenland 
continental  shelf  is  cut  by  three  troughs  or  depressions: 
•Belgica  Trough: 

Belgica  Trough  cuts  the  shelf  from  the  shelf  break  at 
77 °N  westward  to  just  north  of  lie  de  France  at  78 ®N. 

The  greatest  depth  is  in  excess  of  500  m  at  the 
shoreward  end. 

•Norske  Trough: 

Norsks  Trough  runs  from  the  west  end  of  Belgica  Trough, 
to  the  west  end  of  Westwind  Trough  and  parallels  the 
coastline . 

•Westwind  Trough: 

Westwind  Trough  is  located  at  the  north  end  of  the 
Greenland  continental  shelf  and  extends  southeasterly  from 
Engolf’s  Fjord.  It  is  somewhat  shallower  than  Belgica 
Trough . 
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Figure  1.9  The  bathymetric  map  of  the  Greenland  continental 
shelf  was  produced  by  contouring  depths  measured  at  CTD 
stations  during  the  1979  Westwind  and  the  1981  and  1984 
Northwind  cruises  (Bourke  et  al.,  1987). 


II. 


BASIC  THEORY 


A.  "C- VECTOR"  CONCEPT 

1.  Geostrophlc  and  Ageostrophlc  Circulation 

The  ocean  currents,  as  mentioned  in  Chapter  I,  are 
not  truely  geostrophically  balanced.  When  we  compute  the 
velocity  in  the  ocean  by  the  geostrophic  method,  we  depict  the 
major  movement  of  the  oceanic  waters  but  the  total  flow  is  not 
represented,  particularly  in  regions  having  strong  fronts  near 
the  coast  (Arnone  et  al.,  1990).  Neither  geostrophic  currents 
nor  the  dynamic  height  provides  any  information  about  the 
three  dimensional  flow  field  near  oceanic  fronts.  In  order  to 
diagnose  the  three  dimensional  flow  field  near  frontal 
regions,  a  new  theory,  the  'C- vector',  is  introduced.  The 
flow  over  most  continental  shelves  contains  two  parts, 
geostrophic  currents  and  ageostrophic  contributions.  The 
geostrophic  currents  are  obtained  from  hydrographic  data  sets 
while  the  ageostrophic  circulation  includes  the  non¬ 
geos  trophic  contribution,  namely  the  flow  driven  by  the 
surface  wind  field.  The  C-Vector  is  a  new  technique  which  can 
calculate  the  three-dimensional  ageostrophic  circulation  by 
incorporating  hydrographic  and  surface  wind  data  sets  (Chu, 
1992)  . 

Assuming  the  Boussinesq  approximation,  the  geostrophic 
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currents  in  the  north -south  and  east -west  directions  are 


computed  from  the  density  data  derived  from  CTD  observations 
by 

f  g  dp  _  g  _  db  ..  . 

dz  po  dy  po  dy  By 


^  .  g  dfi  _  g  _Bb 

dz  Pg  Bx  Po  Bx  Bx 


(lb) 


where 

Po  :  the  characteristic  value  (a  constant) 
representing  the  mean  density 
p  :  the  density  deviation  from  po  (i.e.,  Po+p  is 
the  observed  density  field) 
a,  :  p  +  Po  -  1000 
b  :  the  buoyancy  force  (*gp/Po) 

As  stated  previously,  the  geostrophic  current  does  not 
represent  the  entire  water  motion.  To  determine  this  we  must 
consider  the  non-geostrophic  contributions.  The  basic  non- 
geostrophic  equations,  describing  the  coastal  water  flow 
without  an  adiabatic  source  of  buoyancy  and  induced  by  the 
wind  stress,  are 

+  u-/(v-yj  (2a) 

ot  ^  dz 

(-^+t^-V)  v+f(u-C7J  (2b) 

ot  ^  dz 


2U 


2j 

~*v-V)  b+N^w=o 
ot 


(2c) 


du  _  dv  _  dw 
dx  dy  dz 


(2d) 


where 


V  :  (u,v,w) 

f  :  Coriolis  parameter  [2Qsin { latitude) ] 

N  :  Brunt -Vaisala  frequency 

(Y*,  Y’')  :  turbulent  momentum  flux  caused  by  the 

surface  wind  stress  (t,,  t^} 


yy| 


(3) 


We  decompose  the  three-dimensional  flow  into  two 
parts:  the  geostrophic  current (Vg)  and  the  ageostrophic 

circulation (VJ 


V  =  V+V^  .  (4) 

V  a 

After  the  decomposition,  the  basic  equations  (2a) - (2d)  should 
have  the  following  form: 

H(i?^)  =  0i(v^,  p,  yy)  (5) 


where  H  and  are  differential  operators.  If  the  real  forms  of 
these  operators  are  obtained,  we  can  use  (5)  to  compute  the 
three-dimensional  ageostrophic  circulation  (V,)  from  the 
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hydrographic  data  (V„,p)  and  the  surface  wind  data. 

2.  The  C- vector 

If  a  quasi -geostrophic  (QG)  system  is  employed,  this 
leads  to 


dt 


V  ■  V 


+  V-V 
dt  ^ 


{6b) 


The  basic  equations  for  the  coastal  water  (2a) - (2d)  on  an  f 
plane  can  be  transformed  to  the  following  equations: 


(7a) 


fu  =  -  ( —  +  i?  •  V)  y 

^  dz  dt  ^  ^ 


{7b) 


=  -  ( -^ + y^-V)  b 


(7c) 


du.  dv.  dw, 

+  +  -rr^  =  0 


dx  dy  dz 


(7d) 


Using  the  cross -derivatives  with  equations  (7a) 
;7c) ,  the  ageostrophic  pseudo-vorticity  equations  become 


(8a) 
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(8b) 


where 


_d_ 

dx 


_d_ 

dy 


=  2C, 

=  2C3 


(8c) 


dU  dV 
C  =  -f{  ^ 

^  ^  ay  dz 


£^_^)  ^  f 
dy  dz  2  dz^ 


(9a) 


^  dz  dx  dz  dx  2  dz^ 


(9b) 


c  f  d  ,dY^  ^  dyy 

^  dx  dy  dx  dy  2  dz^  dx  dy 


(9c) 


Cs  (q  ,  Cj ,  C3)  (9d) 

C, ,  C2  and  C3  are  the  three  components  of  the  C-Vector.  The  C- 
vector  can  be  computed  from  the  oceanic  and  wind  geostrophic 

1 

currents  (Ug,  V^)  which  are  obtained  from  the  hydrographic  data 
set  and  the  wind  data  {Y\  Y’')  . 

In  addition,  we  have  another  important  vector,  the 
pseudo-vorticity  of  total  flow  (C,,  q,  q)  ,  which  is  induced 
from  the  geostrophic  and  the  ageostrophic  currents.  These 
three  components  are  found  by  following  process: 
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We  decompose  v  and  w  into  v  =  +  v,  ,  w  =  Wg  +  but  Wg  =  0 , 

so  w  =  Wj  ,  then  we  get 


c 

*  f2  2  dz 


(10a) 


After  a  similar  operation,  we  can  get  the  other  components,  Cy 
and  ; 


c  3-£k  +2 ^ 
^  f2  2  dy 


(lOb) 


and  the  ageostrophic  pseudo-vorticity  is  most  clearly  seen  in 
the  nondimens ional  form  for  flows  with  constant  f  and  N: 
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VxV^=2R^C  (11) 

V-V^  =  0  (12) 

where  Rg  s  U/(FL)  is  the  Rossby  nximber,  H/L  =  W/U  =  f/N  is 
chosen  for  the  scaling,  and  C  is  nondimensional  here.  The 
ageostrophic  vorticity  is  proportional  to  the  C-vector; 
therefore,  a  C-vector  streamline  can  be  viewed  as  an 
ageostrophic  vortex  line  (Figure  2.1). 


Figure  2.1  Ageostrophic  circulation  induced  by  C  in 
three-dimensional  space. 


b.  Inducing  vertical  velocity  from  the  vertical 


vorticity  of  the  C- vector 

The  vertical  velocity  equation  is  obtained  from 
■k*Vx [eq(12) ]  : 


=-2i?,ic*VxC  (13) 

a  V 

This  indicates  that  the  vertical  velocity  is  induced  by  the 
vertical  vorticity  of  C  (Figure  2.2). 

If  we  define 


Figure  2.2  Vertical  velocity  w,  induced  by  horizontal 
rotation (vorticity)  of  C. 
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then  equation  (14)  becomes 


V^w^=-2R^4>  .  {14b) 

If  the  normal  modes  are  considered,  the  vertical 
velocity  w,  has  the  same  sign  as  Sk. 

c.  Non -divergence  of  the  C- vector 

To  add  the  differential  of  equation  (9a)  with 


respect  to  x,  (9b)  with  respect  to  y  and  (9c)  with  respect  to 
z,  we  can  get  the  following  equation; 


inferred  from  the  vertical  component  of  the  C- Vector. 
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V'C  =  0 


(15) 


which  indicates  that  the  C-vector  is  non-divergent ,  i.e.,  the 
horizontal  divergence  (convergence)  should  be  balanced  by  the 
vertical  convergence  (divergence).  Figure  (2.3)  shows  tnat  the 
vertical  vorticity  of  ageostrophic  flow  can  be  inferred  from 
the  vertical  component  of  the  C-vector. 

B .  WEIGHTING  FUNCTION 

To  compute  the  C-Vector  three-dimensionally  in  the  ocean, 
the  hydrographic  data  on  a  gridded  surface  are  necessary. 
Since  the  CTD  stations  acquired  during  the  Northwind  1984 
cruise  were  not  on  grid  points,  a  weighting  function  is  used 
to  distribute  the  original  CTD  data  onto  the  grid  points.  The 
data,  after  processing  by  this  method,  will  still  retain  the 
original  characteristics  of  this  particular  special  region  of 
the  ocean. 

Barnes  (1964)  used  a  Gaussian  weighting  function  to  develop 
an  objective  analysis  and  filtering  technique  for  application 
to  a  two-dimensional  meteorological  field.  He  modified  this 
method  in  i.973  for  scale  separation.  Maddox (1980)  developed 
this  analysis  technique  to  study  atmospheric  circulations 
associated  with  a  large  complex  of  thunderstorms.  This  scale 
separation  has  been  used  successfully  to  isolate  different 
subsynoptic  meteorological  phenomena.  Tintore  et  al .  (1990) 
used  this  technique  to  investigate  mesoscale  dynamics  and 
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vertical  motion  in  the  Alboran  Sea.  This  technique  is  also 
useful  in  the  analysis  of  oceanographic  phenomena. 

The  two-dimensional  values  of  a  parameter,  e.  g., 
temperature,  f'^(i,j)  at  each  incremental  depth  at  the  grid 
point  (i,j)  are  computed  from  the  observed  data  from  a  total 
of  N  stations.  There  are  four  basic  steps. 

The  first  step  is  to  calculate  the  first  guess  value  at 
each  grid  point  from  the  observed  data  at  each  station. 

.V 

52  j) 

-  (16) 

nml 

where 

fo(x,y)  :  the  observed  data 

f°(i,j)  :  the  first  guess  value  at  each  grid 

point  from  the  observed  data  at  each 
station. 

The  weighting  function  w„(i,j)  is 

W,  (i,  j)  =exp[-d„(i,  j)  V4c]  (17) 

where 

d„(i,j)  :  the  separation  distance  between  the  grid 

point  (i,j)  and  the  nth  station, 
c  :  1000  (analysis  parameter) . 

The  second  step  is  to  calculate  the  first  guess  value  at 
each  station  considering  all  the  stations.  The  weighting 
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function  in  eq(2l)  is 


,V 

{x,y)fjx.y} 

- - -  (18) 

.V 

53  K’  (x-y) 


W„’  (x,y)  =exp[-d„(x,y)^/Ac]  (19) 


where 

d„(x,y)  :  the  distance  between  a  station  and 

all  other  stations. 

The  third  step  is  to  calculate  the  difference  between  the 
original  station  data  and  the  first  guess  value  at  each 
station. 

Af„(x,y)=f„(x,y)-f„°(x,y)  (20) 


where 

fa(x,y) 

C(x,y) 


Af„(x,y) 


the  observed  data  at  each  station, 
first  guess  value  of  tae  nth  station 
determined  from  adjusting  the 
observed  data  at  all  stations, 
the  difference  between  the  observed 
data  and  the  first  guess  value  at 
the  same  point. 
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The  fourth  step  is  to  compute  the  resulting  grid  point 
data  as  following  : 

.V 

(i,  j)Af„{x,y) 

- - -  (21) 

where  the  modified  weighting  function,  w„*’,  is 

j)  =exp[-d„(i,  j)  V4cg]  ,  0<g<l  (22) 

Advantages  of  Using  Weighting  Functions: 

Barnes (1973)  listed  the  advantages  of  this  Gaussian 
weighting,  objective  analysis  and  filtering  technique  as  the 
following  : 

•  The  weighting  function  constants  c  and  g  are  chosen  prior 
to  the  analysis  so  that  pattern  scales,  resolvable  by  the 
data  distribution,  will  be  revealed  to  a  known  response 
amplitude. 

•  If  Wn  approaches  zero  asymptotically,  the  influence  of 
data  may  be  extended  to  any  distance  without  changing  the 
weighting  function  and  response  characteristics. 

•  Small-scale  noise  is  suppressed  so  that  further  smoothing 
with  numerical  filters  is  not  necessary. 

•  Desired  scale  resolution  is  r.chieved  with  only  one 
iteration. 


31 


III. 


DATA  PROCESSING 


As  previously  mentioned,  there  were  333  CTD  stations 
acquired  during  the  1984  Northwind  cruise.  Of  these,  the 
density  (a,)  data  from  237  stations  were  chosen  for  use  in  this 
thesis.  These  data  can  be  used  to  compute  the  cr,  at  each  grid 
point.  After  determining  the  data  at  each  grid  point,  the 
geostrophic  velocity  and  the  C-vector  at  each  grid  point  can 
then  be  computed. 

A.  COMPUTING  THE  ASSUMED  a,  VALUE  AT  EACH  GRID  POINT 

The  weighting  function  described  in  the  previous  chapter 
can  be  used  to  compute  the  assumed  a,  value  at  each  grid  point 
from  the  a,  data  of  the  original  CTD  station.  The  following 
steps  constitute  the  procedures  for  obtaining  the  c,  data  at 
each  grid  point: 

(1)  choose  the  grid  area; 

(2)  set  the  relative  coordinate  for  each  grid  point  and 
each  CTD  station; 

(3)  compute  the  distance  between  each  grid  point  and  each 
CTD  station,  as  well  as  between  each  CTD  station  and  every 
other  CTD  station  from  the  relative  coordinate; 

(4)  obtain  the  'grid  to  station'  weighting  value  from  the 
'grid  to  station'  distance  data  file; 
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(5)  obtain  the  'station  to  station'  weighting  value  from  the 
'station  to  station'  distance  data  file; 

(6)  calculate  the  first  guess  of  the  a,  value  at  each 

grid  point  in  every  layer  by  using  the  'grid  to  station' 
weighting  value; 

(7)  calculate  the  guessed  station  value  from  the 
original  cr,  value  of  the  other  stations  by  using  the 
'station  to  station'  weighting  value; 

(8)  obtain  the  differences  by  subtracting  the  guessed  value 
at  each  CTD  station  from  the  guessed  station  a,  value; 

(9)  obtain  the  assumed  differences  between  the  first 
guessed  value  at  each  grid  point  and  the  resulting 
guessed  cr,  value  by  using  the  weighting  function  to 
calculate  the  differences  in  (8) ; 

(10)  obtain  the  resulting  guessed  ff,  at  each  grid  point  by 
adding  the  first  guessed  <t,  at  each  grid  point  and  the 
assumed  differences  in  (9)  . 

B.  CC3MPUTATION  OP  GEOSTROPHIC  VELOCITY 

The  baroclinic  velocity  can  be  calculated  from  ct,  data  by 
using  equations  (la)  and  (lb)  .  The  constants  in  these 
equation  are  g  =  9.81  m  s'^,  p  o  =  1026  kg  m'^,  the  Coriolis 
parameter  f  =  2nsin0  0  =  7.29xl0‘*  sec*  and  0  =  75°.  The 
geostrophic  method  requires  the  assumption  of  a  level  of  no 
motion.  Because  the  bottom  depths  vary  considerably  over  the 
study  area,  the  weighting  function  can  produce  an  assumed  a, 
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value  below  the  bottom,  that  is,  a  shallow  station  can  refer 
to  the  (T,  value  of  a  deeper  station  around  it  to  assume  the  a. 
value  below  the  bottom  by  calculating  the  weighting  function. 
The  assumed  cr,  value  above  300  m  at  each  station  can  be 
computed.  The  geostrophic  velocities,  Ug  and  Vg,  can  be 
calculated  by  using  the  assumed  cr,  value  at  3  00  m,  chosen  to 
be  the  level  of  no  motion. 


C.  COMPUTATION  OP  THE  C -VECTOR 

After  Ug  and  Vg  are  obtained  by  using  the  geostrophic 
method,  C,,  Cj  and  C3,  the  three  components  of  the  C- vector, 
can  be  computed  continuously  by  using  equations  (9a)  ,  (9b)  and 
(9c)  .  The  turbulent  momentum  fluxes  Y*  and  are  calculated 
from  the  following  equations: 


yx 


where  the  wind  stress  t*  and  is  calculated  from  : 


where  the  drag  coefficient  Co  =  1.3  x  lO'^  and  the  air  density 
p„=1.293  X  10'^  kg  m'^.  |V»1  is  the  average  wind  speed  at  height 
10  m,  Uw  and  v*  are  the  components  of  wind  speed  in  the  x  and 
y  directions  at  each  CTD  station.  Using  the  weighting  function 
as  before,  u^  and  v*  at  each  grid  point  can  be  computed.  The 
average  depth  of  the  mixed  layer  in  the  East  Greenland  Current 
during  the  cruise  of  Northwind  1984  is  about  15  m.  We  assume 
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that  the  turbulent  stress  linearly  decreases  with  depth  in  the 
mixed  layer,  and  there  is  no  stress  below  the  mixed  layer. 
Under  these  conditions,  the  C- vector  with  wind  can  be 
computed.  After  obtaining  the  C, ,  Cj,  Cj  components,  the  total 
C-Vector  (C,,  Cy,  C^}  can  be  calculated  continuously  by  using 
equations  (lOa)  ,  (10b)  and  (10c)  .  Once  C,,  Cy  and  C,  have  been 
obtained,  the  stream  function  (Sk)  can  be  computed  easily  by 
using  equation  (11) . 
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IV.  RESULTS 


A.  THE  COORDINATE  SYSTEM 

The  coordinate  system  is  chosen  such  that  the  origin  is 
located  at  (76°3Q’N,  16°30’W);  the  X-axis  is  zonal  and 

positive  eastward;  the  Y-axis  is  meridional  and  positive 
northward;  and  the  Z-axis  is  vertical  and  positive  upward. 
There  are  10  sections  along  the  X-axis (Sections  A  to  J) ,  21 
sections  along  Y-axis (Sections  1  to  21)  and  31  layers  in  the 
vertical  (each  layer  is  10  m  thick).  Figure  1,8  shows  this 
coordinate  system, 

B.  C- VECTOR  PLOTS 

1.  C(/f^  and  Qjt  In  the  X  and  Y  direction 

The  solid  contours  of  Cjt  and  C,/f^  are  the  zonal 
positive  total  and  ageos trophic  nondimens ional  pseudo 
vorticities,  which  means  the  rotational  axis  is  towards  the 
east.  The  dash  contours  of  C,/f  and  C,/f^  are  the  negative 
total  and  ageostrophic  nondimens ional  pseudo  vorticities, 
which  indicates  the  rotational  axis  is  towards  the  west. 
Appendix  I  shows  these  plots.  Similar  to  the  above,  contours 
of  Cy/f  and  Cj/f^  represent  longitudinal -directed  vorticities, 
i.  e.,  rotational  axies  are  pointed  towards  the  north  (solid 
lines)  or  south  (dashed  lines) .  These  are  shown  in  Appendix 
II. 
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2 .  ^  plots 

The  vorticity  of  the  C-vector  is  plotted  on  the 
horizontal  plane  for  each  10  m  layer  in  the  study  area.  The 
function  has  the  same  sign  as  the  ageostrophic  vertical 
velocity.  Therefore,  the  positive  contours  (solid  line)  mean 
upwelling  and  the  negative  contours  (dash  line)  mean 
downwelling.  Appendix  III  shows  these  plots. 

C.  SIGNIFICANT  FEATURES  FOUND  FROM  THE  C-VECTOR  METHOD 

1.  Cross -Coastal  Circulation  Generated  by  the  Surface 

Wind 

The  wind  data  collected  from  the  1984  Northwind  cruise 
(Figure  1.7)  indicates  the  average  wind  speed  is  5  m/s  from 
NNE  to  NE.  The  Ekman  transport  is  directed  to  the  right  of  the 
wind.  From  the  continuity  equation  such  a  wind  field  will 
generate  coastal  downwelling.  This  indicates  that  Cj/f^  should 
be  negative  near  the  coast.  Figure  4.1  shows  the  relation 
between  the  wind  direction  and  the  vertical  circulation  near 
the  east  Greenland  coast. 

Figure  4.2  shows  three  sections  (7,  8  and  9)  of  Cj/f^ 
contours  in  three  dimensional  view.  The  east  Greenland  coast 
lies  to  the  negative  horizontal  vortices  are  found  to  the  west 
in  these  sections,  indicating  that  downwelling  was  occurring 
near  the  coast  and  upwelling  was  occurring  over  the 
outershelf,  strongest  near  the  shelf  break. 
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Figure  4.1  The  relation  between  the  wind  direction  and  the 
vertical  circulation  near  the  east  Greenland  coast. 

From  the  plots  of  C^/f^  in  Appendix  II,  we  can  see  that 
horizontal  rows  appear  from  the  surface  to  40  m  depth.  This 
indicates  that  wind  effects  dominate  the  vertical  circulation 
in  the  upper  ocean. 

2 .  Vertical  Circulation  Induced  by  the  Antlcyclonlc  Gyre 

A  northward  flow  in  Norske  Trough  results  from  a 
clockwise  turning  of  the  current  around  the  southwest  corner 
of  Belgica  Bank  (Tunniclif fe,  1985) .  Figure  4.3  shows  the 
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relative  position  between  the  grid  points  and  the  surface 
circulation  in  EGC.  The  anticyclonic  gyre  is  found  at  the 
center  the  grid  area.  The  vertical  motion  should  be  upward  at 
the  gyre  center,  because  the  water  in  the  center  of  gyre  will 
move  out  from  the  gyre  center  to  conserve  continuity.  Hence 
upwelling  should  occur  in  the  central  part  of  the  gyre.  Figure 
4.4  shows  three  dimensional  Cy/f  contours  for  Sections  10,  ll 
and  12. The  dashed  contours  indicate  cyclonic  or  upward  motion; 
conversely  for  the  solid  contours.  The  central  part  of  the 
gyre  is  located  from  X  =  75  km  to  X  =  125  km  in  section  10, 
from  X  =  50  km  to  X  =  175  km  in  section  11,  from  X  =  75  km  to 
X  =  175  km  in  section  12.  We  can  see  the  surface  upwelling  in 
this  region,  which  suggests  that  the  upwelling  is  induced  by 
the  presence  of  the  anticyclonic  gyre. 

The  upwelling  induced  by  the  anticyclonic  gyre  is  also 
seen  in  the  field.  The  central  part  of  this  anticyclonic 
gyre  is  located  at  X  =  50  to  150  km  and  Y  =  200  to  325  (Figure 
4.5) .  The  ^-values  at  z  =  0,  10,  and  20  m  in  the  central  part 
of  the  gyre  are  generally  positive  (upward  motion) ,  which  also 
confirms  the  results  obtained  from  Cy- field. 

3 .  Vertical  Circulation  Generated  by  the  EGPF 

The  EGPF  was  investigated  in  previous  studies.  Figures 
1.3  and  1.4  (Bourke,  et  al . ,  1985)  reveal  the  location  and 
structure  of  the  EGPF  from  the  temperature  and  salinity 
contours.  From  these  two  figures  we  can  see  that  a  sharp 


41 


Figure  4.4  Contours  of  C^/f  at  Sections  10,  11,  and  12.  The 
dashed  line  indicates  upward  motion  and  the  solid  line 
downward  motion. 


temperature  and  salinity  gradient  is  present  below  40  m  (the 
effect  region  by  wind)  and  extends  from  about  60  m  to  120  m. 

In  Figure  4.3  a  strong  frontal  jet  flows  through 
Sections  1  and  2  in  the  Y  direction.  C^/f-*  contour  plots  for 
these  two  sections,  shows  that  a  two* cell  structure  is  found 
at  a  depth  about  50  to  120  m  near  the  EGPF.  Therefore,  this 
two- cell  structure  indicates  that  the  AIW  of  the  RAC  meets 
with  PW,  a  feature  which  can  be  seen  in  X-Z  cross  sections  1 
and  2  at  a  depth  of  about  60  to  120  m  {Figure  4.6) .  Figure 
4.7  shows  that  a  two- cell  structure  occurs  at  the  front 
between  AIW  (warm  and  saline)  and  the  PW  (cold  and  fresh) .  We 
can  sketch  a  frontal  model  for  the  EGPF.  When  AIW  and  PW  are 
mo  ing  Coward  each  other,  two  horizontal  opposite  vortices 
will  be  formed  with  a  downward  motion  along  the  front. 
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Figure  4.6  Three  dimensional  Cj/f^  contours  at  Section  1,  2, 
3.  The  EPGF  passes  through  Sections  1  and  2. 


Figure  4.7  A  two- cell  structure  with  downward  motion  occurs 
in  the  north- south  direction,  AIW  {warm  and  saline;  at  the 
right,  PW  (cold  and  fresh)  at  the  left. 
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V.  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  SUMMARY 

A  new  technique,  the  'C-vector  method',  is  employed  to 
define  the  three-dimensional  circulation  of  the  East  Greenland 
Current  from  CTD  and  wind  data  collected  during  the  Northwind 
1984  cruise  to  the  Greenland  Sea. 

Three  different  types  of  vertical  circulation  can  be 
identified  by  the  C-vector  method;  (1)  the  wind-driven  coastal 
circulation,  (2)  an  anticyclonic  gyre  which  induces  vertical 
circulation,  and  (3)  a  front  induced  circulation.  The  results 
demonstrate  that  the  C-vector  method  is  useful  in  calculating 
the  ageostrophic  circulation  using  CTD  and  wind  data.  The 
ageostrophic  flow  can  sometimes  be  significant  and  the  purpose 
of  this  study  is  to  particularly  emphasize  the  importance  of 
the  ageostrophic  circulation  in  the  waters  off  the  east 
Greenland  coast  as  inferred  by  the  C-vector  using  CTD  and  wind 
data . 

B.  CONCLUSIONS 

•  The  C-vector  method  is  a  very  good  tool  to  diagnose  the 
three-dimensional  pseudo-vorticity  fields  in  the  ocean  by 
using  CTD  measurements  and  wind  stress  data. 

•  The  vertical  motion  in  the  ocean  can  be  investigai *  J  from 
the  vertical  vorticity  of  C-vector  for  the  total  flow,  value. 
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•  If  the  quasi -geostrophic  assumption  can  be  removed  in 
future  work,  more  realistic  results  can  be  obtained. 

C.  RECOMMENDATIONS 

(1)  A  shortcoming  in  the  calculation  of  the  C- vector  was 
the  uncertainty  associated  with  the  determination  of  a 
reference  level  or  level  of  no  motion.  The  determination  of 
a  reference  level  still  remains  disputable  because  there  is  no 
universal  method  which  guarantees  this.  Density  differences 
between  stations  in  the  Greenland  Sea  study  area  demonstrate 
that  300  m  is  a  representative  level  of  no  motion  for  the  deep 
waters.  While  the  density  difference  stations  is  not  always 
zero  at  300  m,  there  is  a  good  indication  that  at  this  level 
the  currents  are  very  small. 

(2)  In  this  study,  we  didn't  solve  the  dynamical 
equations.  In  the  future,  the  pseudo-vorticity  equations  (8a) 
to  (lOc)  should  be  solved  under  certain  boundary  conditions  to 
obtain  the  three  dimensional  ageostrophic  velocity  field. 
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APPENDIX  I  ;  THE  PSEUDO  VORTICITY  CONTOUR  PLOTS  IN  X 
DIRECTION 
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Figure  1.1  C,/f^  contour  at  Section  A  and  section  B 


Figure  1.3  contour  at  Section  E  and  Section  F. 
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Figure  1.9  C,/f  contour  at  Section  G  and  Section  H. 
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Figure  1. 10  C,/f  contour  at  Section  I  and  Section  J. 
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APPENDIX  II  :  THE  PSEUDO  VORTICITY  CONTOUR  PLOTS  IN  Y 
DIRECTION 
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Figure  II. 1  Cj/f^  contour  at  Section  1  and  Section  2 
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Figure  II. 3  Cj/f^  contour  at  Section  5  and  Section  6. 
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Figure  II. 4  Cj/f^  contour  at 
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Figure  11.12  Cy/f  contour  at  Section  1  and  Section  2. 
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Figure  11.13  C^/f  contour  at  Section  3  and  Section  4 
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Figure  11.15  C^/f  contour  at  Section  7  and  Section  8 
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Figure  11.16  C^/f 
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Figure  11.20  Cy/f  contour  at  Section  17  and  Section  18. 
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Figure  11.21  Cy/f  contour  at  Section  19  and  Section  20 
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Figure  11.22  C^/f  contour  at  Section  21. 


82 


APPENDIX  III  i  'If  CONTOUR  PLOTS 


83 


8 


OlWCCTIONlKHi 


85 


LIST  OP  REFERENCES 


Aagaard.  K.,  and  L.  K.  Coachman,  -he  Ease  Greenland  Current 
north  of  the  Cenmark;  Strait,  I,  Arctic,  21  3}  ,  131-2CC,  l?68a. 

Aagaard.  K.  ,  and  L.  K.  Caachanan,  The  East  Greenland  lurrent 
north  of  the  Denmark  Strait,  II,  Arctic,  21 '4!  ,  2c7-23C, 
1968b. 

Arnone,  R.  A.,  D.  A.  Wiesenburg,  and  K.  D.  Saunders,  The 
Origin  and  Characteristics  of  the  Algerian  Current,  J. 
Geophys .  Res.,  95,  1587-1593.,  1990. 

Barnes,  S.  L.,  A  Technique  for  Maximizing  Details  in  -h-a.mericai 
Weather  Map  Analyses,  J.  Appl .  Meteor.,  3.  396-409.,  1964. 

Barnes.  S.  L.,  Mesoscale  Objective  Map  Analysis  Using  Weighted 
Time  Series  Observations,  NOAA  Tech.  Memo.  ERL  NSSL-52,  60  pp . 
[NTIS  COM- 73- 10781] ,  1973. 

Bourke,  R.  H.  ,  J.  L.  Newton,  R.  G.  Paquette,  and  M.  D. 
Tunnicliffe,  Circulation  and  Water  Masses  of  the  East 
Greenland  Shelf,  J.  Geophys.  Res.,  92(C7),  6729-6740,  1987. 

Chu,  P.,  C.,  Three  Dimensional  Eastern  Greenland  Sea 
Circulation  Compute  from  a  CTD  Data  Set,  ARCSS  Ocean- 
Atmosphere-  Ice  Interactors,  NSF,  pp .  61-64,  1992. 

Chu,  P.,  C.  ,  Three  Dimensional  Pseudo- Vorticity  Fields  in  the 
West  Spitsbergen  Current,  Polar  Meteorology  and  Oceanograohy , 
vol.  Ill,  pp.  101-104,  1992. 

Hoskins,  B.,  J.,  Draghich  I.  and  Davies  H. ,  C.,  A  New  Look  at 
the  Q-equation,  Quart.  J.  Roy.  Meteor.  Soc.,  104,  31-33,  1978. 

Johnson,  G.  L.,  and  0.  B.  Eckhoff,  Bathymetry  of  t.he  North 
Greenland  Sea,  Deep-Sea  Res.,  13,  1161-1173,  1966. 

Maddox,  R.  A.,  An  Objective  Technique  for  Separation 
Macroscale  and  Mesoscale  Features  in  Meteorological  Data, 
American  Meteorological  Society,  108,  1108-1121,  1980. 

Muench,  R.  D.,  G.  S.  E.,  Largerloef ,  and  J.  T.  Gunn,  1984-85 
Current  Observations  in  the  East  Greenland  Current:  A 
Preliminary  Description,  MIZEXBull.  7,  CRREL  Spec.  Rep.  86-3, 
pp.  41-53,  U.S.  Army  Cold  Reg.  Res.  and  Eng.  Lab.,  Hanover,  N. 
H.,  1986. 


86 


Paquette,  R.  G.,  R.  H.  Bourke,  J.  L.  Newton  and  w.  ?.  Perdue, 
The  East  Greenland  Polar  Front  in  Autumn,  J.  Geophys .  Res., 
90{C3),  4866-4882,  1985. 

Perry,  R.  K.,  H.  S.  Flemming,  N.  2.  Cherkis,  R.  H.  Feden,  and 
P.  R.  Vogt,  Bathymetry  of  the  Norwegian-Greenland  and  Western 
Barents  Seas,  U.  S.  Naval  Research  Laboratory  -  Acoustics 
Division,  Environmental  Sciences  Group,  William  and  Geintz  Map 
Corporation,  Washington,  D.C.  1980. 

Sutcliffe,  R.  ,  C.  ,  A  Contribution  to  the  Probie.m  of 
Development,  Quart.  J.  Royal.,  Meteor.  Soc.  73,  pp  370-383, 
1947. 

Tintore,  J. ,  D.  Gomis,  S,  Alonso,  and  G.  Parrilla,  Mesoscale 
Dynamics  and  Vertical  Motion  in  the  Alboran  Sea,  American 
Meteorological  Society,  811-823,  1990. 

Tunniclif fe,  M.  K.,  An  Investigation  of  the  Waters  of  the  East 
Greenland  Current,  Master's  thesis,  136  pp.,  Naval  Post  grad. 
Sch.,  Monterey,  Calif.,  1985, 

Xu,  Q.,  Ageostrophic  Pseudovorticity  and  Geostrophic  C-Vector 
Forcing  a  New  Look  at  the  Q-Vector  in  Three  Dimensions, 
Journal  of  Atmospheric  Science,  49,  981-990,  1992. 


87 


INITIAL  DISTRIBUTION  LIST 


Center,  Applied  Ahystcs  Labcrat 
Was hing t cn  r  3 1 C  5 

2.  Curry,  Judith 

Department  of  Aerospace  Engineering  Science 
University  of  Colorado 
Campus  Box  429 
Boulder,  Colorado  30305 


Colony,  Roger 
Polar  Science 
University  of 


5.  Office  of  Naval  Research  ; Code  1125  ATJ 
Environmental  Science  Directorate 
800  N  Quincy  Street,  Code  1125 
Arlington,  Virginia  22217 

4.  Defense  Technical  Information  Center 
Cameron  Station 

Alexandria,  VA  22304-5145 

5.  Library,  Code  52 

Naval  Postgraduate  School 
Monterey,  CA  93943-5002 

6.  LCDR.  Wei-Szu  Li 

21-6f  Ruey-Fong  St.,  Ku-San  Dist.,  Kaohsiung, 
Taiwan,  R.O.C. 

7 .  Library 

Chung  Cheng  Institute  of  Technology 
Tau-Yuen,  Taiwan,  R.O.C. 

8 .  Department  of  Oceanography 
Naval  Postgraduate  School 
Monterey,  CA  93943 -.5000 
Attn:  Chairman,  (Code  OC/Co) 

Attn:  Dr.  P.  C.  Chu  (Code  OC/Cu) 

9.  Commanding  Officer 
Navy  Headquarters 
Taipei,  Taiwan,  R.O.C. 

10.  Commanding  Officer 

Naval  Hydrographic  &  Oceanographic  Office 
Tso-Ying,  Kaosh.Lung,  Taiwan,  R.O.C. 


88 


11.  Commanding  Officer 

Naval  Oceanographic  Iff 

NSTL  Station 

Bay  St.  Louis,  MS  39b22 


12.  Prof.  Untersteiner ,  Norherc 
Department  of  .^.cmcspheric  Sciences  .PK--40 
University  of  Washington 

Seattle,  Washington  ?319E 

13.  Prof.  Walsh,  John  E. 

University  of  Illinois 
Department  of  Atmospheric  Sciences 
105  S.  Gregory  Avenue 

Urbana,  Illinois  6l8Qi 


89 


